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LOW-PRESSURE TURBINE SEPARATION CONTROL - COMPARISON WITH EXPERIMENTAL DATA 
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ABSTRACT 

The present work details a computational study, using the Glenn- 
HT code, that analyzes the use of vortex generator jets (VGJs) to 
control separation on a low-pressure turbine (LPT) blade at low 
Reynolds numbers. The computational results are also compared with 
the experimental data of Bons et al. [1] for steady VGJs. It is found 
that the code determines the proper location of the separation point on 
the suction surface of the baseline blade (without any VGJ) for 
Reynolds numbers of 50,000 or less. Also, the code finds that the 
separated region on the suction surface of the blade vanishes with the 
use of VGJs. However, the separated region and the wake 
characteristics are not well predicted. The wake width is generally 
over-predicted while the wake depth is under-predicted. 

NOMENCLATURE 

B jet blowing ratio = (pu) je /(pu) local 

C p blade pressure coefficient = (p oin - p s )/0.5(pu 2 ) in 

C x axial chord of the blade 

d injection hole diameter 

k turbulence kinetic energy 

Q turbulence length scale 

p pressure 

Pr Prandtl number 

Re inlet Reynolds number = C x (pu/p) in 

Tu freestream turbulence intensity 

u streamwise mean velocity 

v* shear velocity 

x,y Cartesian coordinate system with origin at the blade leading 
edge 

y + distance in wall coordinates = pyv*/p 

a thermal diffusivity 

Ay distance (from the wall) of the first point off the wall 
e turbulence dissipation rate 

y wake loss coefficient = (p o in - p 0 , ex )/(p 0iin - p s ,J 
p dynamic viscosity 

p density 

T shear stress 


CO specific turbulence dissipation rate = e/k 

Subscripts 

ef effective value 

ex value at cascade exit 

in value at cascade inlet 

jet vortex generator jet conditions 

C laminar value 

local local blade mid-channel conditions 
o stagnation value 

s static value 

t turbulent value 

w value at the wall 

INTRODUCTION 

For the flow over an airfoil, there is an adverse pressure gradient 
region over a large part of the suction surface. This can cause the flow 
to separate on the suction surface, especially at low Reynolds numbers, 
in which case the boundary layer over the blade may be laminar. It is 
well known that a laminar boundary layer can withstand only a minor 
adverse pressure gradient before separation while a turbulent boundary 
layer can withstand a much stronger adverse pressure gradient region. 
Flow separation reduces lift and increases drag. It is thus undesirable 
and needs to be controlled. 

Flow separation can be controlled by both passive and active 
techniques. The former employs small vortex generators (rectangular 
or delta shaped winglets) which are imbedded in a boundary layer 
ahead of a line of flow separation [2], Spanwise arrays of vortex 
generators are often placed along a wing upstream of the flap hinge or 
inside the lip of a jet engine inlet diffuser. The longitudinal vortices 
generated in the boundary layer increase cross-stream mixing of 
streamwise momentum and suppress or eliminate separation. These 
fixed vortex generators have the advantages of simplicity, ruggedness 
and low cost. However, they add parasitic drag in flow situations 
where separation control is not needed, e.g., take-off and landing. 


*NAS A Resident Research Associate at Glenn Research Center. 


NASA/CR— 2002-21 1689 


1 


In one of the active techniques, small jets blown through holes in 
the solid surface can generate longitudinal streamwise vortices in a 
boundary layer. These vortices increase cross-stream mixing of 
streamwise momentum. The surface holes are pitched an angle to the 
surface (generally low pitch angle of 30-45°), and inclined relative to 
the main flow direction (skew angle varying from 45-90°), much like 
the compound-angled holes on the shower-head of film-cooled blades. 
This method is called the vortex generator jet (VGJ) method. This 
method was first examined by Wallis [3, 4], and by Wallis and Stuart 
[5] primarily for the purpose of delaying shock-induced separation of 
turbulent boundary layers. Ball [6] employed VGJs alone and together 
with fixed generators for stall suppression in inlet diffusers of jet 
engines. Henry and Pearcey [7] show clearly the development of a 
single, dominant vortex for a VGJ at high skew angles. It has been 
shown by Johnston and Nishi [8] that this vortex energizes the 
separating boundary layer by bringing high momentum freestream fluid 
down to the wall. Experimental data on the beneficial effects of VGJs 
are available in [9]. 

The flow requirements for the jets in the VGJ method is a very 
small fraction of the total jet engine flow, and compressor bleed air is 
sufficient to power the system. The jets can be turned on and off as 
desired. With jets off, there is no parasitic drag. The potential 
applications of the VGJ method are numerous for external flow over 
aircraft and missiles. Internal flows in diffusers can also be improved 
and pressure losses reduced. The VGJ method can be used with 
appropriate flow direction sensors and a feedback control system to 
suppress surge and rotating stall in compressors of jet engines. 

The operating Reynolds number for the low-pressure turbine in an 
aircraft gas turbine engine can drop below 25,000 during high altitude 
cruise. At such low Reynolds numbers, the boundary layers on the 
LPT blades are largely laminar, making them highly susceptible to 
flow separation near the aft portion of the blade suction surface, with 
associated increase in losses and drop in performance. Altering the 
blade shape to avoid this low Reynolds number separation problem is 
not feasible since such a modification will likely impair the engine 
operation at higher (design) Reynolds numbers. As such, an active 
flow control strategy, such as the use of VGJs, has been tested 
experimentally at the Air Force Research Laboratory recently [1, 10], 
This paper aims to use the Glenn-HT code in order to analyze the use 
of VGJs at low Reynolds numbers, and to compare the computational 
predictions with the experimental data of Bons et al. [1], 

ANALYSIS 

The numerical simulation has been performed using the NASA 
Glenn Research Center General Multi-Block Navier-Stokes Convective 
Heat Transfer code, Glenn-HT. Briefly, the code, formerly known as 
TRAF3D.MB [11], is an explicit, multigrid, cell-centered, finite 
volume code with a k -(0 turbulence model without any wall functions. 
This is a general purpose flow solver designed for simulations of flows 
in complicated geometries. The Navier-Stokes equations in a rotating 
Cartesian coordinate system are mapped onto a general body-fitted 
coordinate system using standard techniques. The multistage Runge- 
Kutta scheme developed by Jameson et al. [12] is used to advance the 
flow solution in time from an initial approximation to the steady state. 
A spatially varying time step along with a CFL number of 4 is used to 
speed convergence to the steady state. Eigenvalue-scaled artificial 
dissipation and variable-coefficient implicit residual smoothing are used 
along with a full-multigrid method. The overall accuracy of the code 
is second order. No wall functions are used, thus avoiding any bias to 


the complex three-dimensional flow structures near the blade or any 
other surface. While the Glenn-HT code is basically a compressible 
code, the present problem involves incompressible flow at rather low 
Reynolds numbers. Thus, there were some questions regarding the use 
of Glenn-HT code in such a situation. Needless to say, such questions 
about getting converged solution from the Glenn-HT code for such a 
case have been put to rest. How these results compare with 
experimental data is another matter. 

While the Glenn-HT code has the original k - to model [13], the 
shear stress transport (SST) model of Menter [14], and the k - to model 
of Wilcox [15] were implemented in it by Garg and Ameri [16]. The 
SST model encompasses both the k -(0 and the k-e models, with the 
original k - co model of Wilcox [13] activated in the near-wall region 
and the standard k-£ model [17] activated in the outer wake region and 
in free shear layers. Moreover, the definition of eddy viscosity is 
modified to account for the transport of the principal turbulent shear 
stress. The reader is referred to Menter [14] for an elucidating 
discussion of the SST model. More details on the relevant equations 
and their implementation are available in Garg and Ameri [16]. The 
SST model was used for the present computations. 

It is assumed that the effective viscosity for turbulent flows can 
be written as 

+ (1) 

where the laminar viscosity p ( is calculated using a power-law for its 
dependence on temperature [18], The turbulent viscosity p t is 
computed using the SST model. The turbulent thermal diffusivity is 
computed from 

a t = - A. (2) 

P Fr , 

where a constant value of 0.9 is used for the turbulent Prandtl number, 
Pr,. 

Boundary Conditions 

At the main flow inlet boundary located at an axial distance equal 
to the blade axial chord upstream of the blade leading edge, the total 
temperature, total pressure, whirl, and meridional flow angle are 
specified, and the upstream-running Riemann invariant based on the 
total absolute velocity is calculated at the first interior point and 
extrapolated to the inlet. The velocity components are then decoupled 
algebraically, and the density is found from total temperature, total 
pressure and total velocity using an isentropic relation. For the 
turbulence model, the value of k and co is specified using the 
experimental conditions, namely 

k = l.5(u in Tu in ) 2 , a>=fc 1/2 /e, (3) 

where Tu m is the intensity of turbulence at the inlet (taken to be 0.01 
as per experimental data for the PakB blade), u in is the absolute 
velocity at inlet, and C is the integral length scale representing the size 
of the energy containing eddies. This length scale needs to be 
assumed, if not reported as part of the experimental conditions, as in 
the present case. It was assumed to be 5% of the blade axial chord. 

At the main flow exit plane located at an axial distance equal to 
80% of the blade axial chord downstream of the blade trailing edge, 
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the static pressure is specified and the density and velocity components 
are extrapolated from the interior. At the solid surface of the blade, 
the no-slip condition is enforced, and temperature is specified as per 
experimental data. The boundary conditions for turbulence quantities 
on the walls are k = 0, and 


w = 100 — 

dy 


vail 


(4) 


for a hydraulically smooth surface. An upper limit is imposed on the 
value of (D at the wall, as suggested by Menter [19] and found 
effective by Chima [20], 


(oj ) 
v max-' w 


800 v 

Re (Ay ) 2 


(5) 


The grid around the blade extends to mid-way between two 
adjacent blades with periodic flow conditions in terms of cylindrical 
velocity components set on a dummy grid line outside this boundary. 
For a linear cascade (which is true for the experimental data), it is 
possible to consider only a slice of the real span for computational 
purposes with a periodic boundary condition at both ends of the 
computational span. For the basic blade with no VGJ, this slice can 
be of any reasonable width, but with VGJs, the width of the slice is 
one span-wise pitch of the holes. 

The effect of VGJs has been incorporated in the form of 
appropriate boundary conditions at the hole locations on the blade 
surface. Each hole exit in its true oval shape is represented by 224 
control volumes. Different velocity and temperature profiles for the 
injected fluid can be specified at the hole exits. For the present study, 
polynomial distribution [21] of velocity and temperature of the 
injectant at the hole exit was specified. The injectant angle was taken 
to be the same as the hole angle. Turbulent intensity at the hole exit 
was assumed to be 1% (same as for the freestream), while the 
turbulence length scale at the hole exit was taken to be the same as the 
hole diameter d. 


EXPERIMENTAL DETAILS 


Measurements were made by Bons et al. [1] in an induction wind 
tunnel at the U.S. Air Force Research Laboratory. The linear cascade 
consisted of eight 0.88 m (34.5 in) span, 0.18 m (7 in) axial chord (C x ) 
blades plus two partial blade endwalls. The blades were fabricated 
from molded polyurethane resin. The 2-D blade shape studied is the 
Pratt & Whitney "PakB" research design, which is a Mach number 
scaled version of a typical highly loaded LPT blade design. The 
cascade had a solidity (axial chord to blade spacing) of 1.13, an inlet 
flow angle of 55° (measured from the plane of the cascade), and a 
design exit angle of 30°. Blades 4 and 6 were instrumented around 
both pressure and suction surfaces (in the center 0.2 m of their span) 
with forty 1 mm diameter static pressure taps. Uncertainties in the 
pressure measurement translated to an error of ± 0.18 in the C p data at 
Re = 100.000 [1], 

Blade #5 of the cascade had the active separation control. It was 
manufactured with a hollow cavity running the span of the blade and 
covering the region from 40% to 90% axial chord. A needle valve 
located upstream of the feed port allows fine control of the mass flow 
rate into the blade cavity. Air exhausts from this valve into a 1.2 cm 
diameter capped copper tube running the span of the blade inside the 
cavity. Holes of 1.5 mm diameter and spaced every 2.54 cm along the 


copper tube produced an even distribution of air flow to the VGJs. 
The 1 mm diameter (d) cylindrical VGJ holes were drilled from the 
exterior surface of the blade, through the 4 mm thick wall, and into the 
cavity with a 30° pitch angle and a 90° skew angle. The VGJ holes 
have a length of 8 d and are spaced every 10c/ along the center 0.46 m 
of this blade span. The row of VGJs was placed at 0.73 C x . This 
location is within the separated region on the baseline, uncontrolled 
blade at a Reynolds number of 50,000 [1]. The jet blowing ratio (B) 
was computed as the ratio of the jet exit velocity to the local 
freestream velocity as calculated from the local pressure coefficient. 
While the inlet Reynolds number was varied between 25,000 and 
100.000, B was taken as 0, 1, 2 or 4. 

Standard pitot and hot-wire probes mounted on a large Dantec 3- 
axis traverse located atop the tunnel facility were used to map a 0.6m 
x 0.6m planar cross-section of the flow domain. Wake velocity 
profiles were measured using a single 4pm diameter tungsten hot-wire. 
A single TSI sub-miniature hot-film probe was used for boundary layer 
profile measurements. The hot- film sensor diameter was 25pm and 
the sensing element length was 0.25mm. The error in the hot-wire and 
film probes was within ±2% at flow rates of interest. More details are 
available in [1], 

COMPUTATIONAL DETAILS 

The computational span extended over a part of the blade span 
with a periodic boundary condition at both ends. In the axial direction, 
the computational domain extended from the inlet plane located one 
axial chord upstream of the blade leading edge to the exit plane located 
80% of the axial chord downstream of the blade trailing edge. Around 
the blade, the grid extends to mid-way between two adjacent blades 
with periodic boundary conditions. Figure 1 shows a spanwise section 
of the multi-block viscous grid around the blade. The viscous grid is 
obtained from an inviscid grid by clustering the grid near all the solid 
walls (blade here). The clustering is done in such a way as to ensure 
that in the viscous grid, the distance of any cell center adjacent to a 
solid wall, measured in wall units (y + ), is less than half for the cases 
studied here, following Boyle and Giel [22]. The inviscid grid was 
generated using the commercial code GridPro/az3000 [23], For 
computational accuracy the ratio of two adjacent grid sizes in any 
direction was kept within 0.8-1.25. As can be observed from Fig. 1, 
the grid quality is very good especially near the blade surface. Figure 
1 also shows the grid over the blade span near and within the injection 
hole exits. 

Initially, the grid for the blade with injection holes consists of 136 
blocks but before the solver is used, it can be merged into just 10 
blocks using the Method of Weakest Descent [24], The final viscous 
grid consists of 160128 cells, formed by clustering near the blade from 
an inviscid grid with 92928 cells. There are 224 cells within the hole 
exit on the blade surface. The inviscid grid has 144 cells around the 
blade (for the O-grid around the blade), 28 cells in the blade-to-blade 
direction from the blade to the periodic boundary in-between the two 
blades, and 24 in the spanwise direction. After clustering, the number 
of cells in the blade-to-blade direction increases to 48. Two more 
grids were generated for a grid-independence study. One inviscid grid 
had 1.5 times the number of cells in each direction as compared to the 
basic grid described above. For the second grid, the basic inviscid grid 
was clustered near the blade with a grid spacing half of that for the 
basic viscous grid. All these variations of the basic grid yielded nearly 
the same values for the skin friction coefficient on the entire blade 
surface as the basic grid; any variations were within ±2%. The results 
presented here correspond to the basic grid shown in Fig. 1. 
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Computations were run on the 16-processor C90 supercomputer 
at NASA Ames Research Center. The code requires about 26 million 
words of storage with all blocks in memory for the case with injection, 
and takes about 7 s per iteration for two levels of multi-grid. A case 
requires about 4000 iterations to converge for a case with injection. 
This unusually large number of iterations for convergence is due to the 
computation of incompressible flow at low Reynolds numbers. 
Recalling that the Glenn-HT code is essentially a compressible code, 
the code was initially run to near-convergence for one case with and 
without injection for an exit Mach number of 0.6. The exit Mach 
number was then subsequently reduced in steps of 0.1 to a final value 
of 0.2, which is essentially an incompressible flow condition. This 
obviated the use of any low Mach number preconditioning. For 
subsequent cases, results for this case served as the starting point. 

RESULTS and DISCUSSION 

Several experimental cases for the PakB blade, available in [1], 
were analyzed for comparison. These cases cover three inlet Reynolds 
number values, 25000, 50000 and 100000 for a freestream turbulent 
intensity of 1%. Most of the data in [1] are for Tu = 1%. Values of 
B for which computations were done are 0, 1 and 2. Figure 2 shows 
the computed shear stress on the basic blade surface, without any VGJ, 
for three values of the inlet Reynolds number. In this figure, the zero 
shear stress line is shown for reference. Clearly, the shear stress is 
negative beyond the separation point on the suction surface of the 
blade. The separation point on the suction surface, corresponding to 
T w = 0, is little affected by the Reynolds number; it is at 0.7 1 C x for Re 
= 50,000 and 100,000, but at 0.73 C x for Re = 25,000. There appears 
to be a small region on the pressure surface as well that has separated 
flow for Reynolds numbers > 50,000. The small separated region on 
the pressure surface is over 0.22 < ,t/C x < 0.25 for Re = 50.000, and 
over 0.18 < .v/C x < 0.29 for Re = 100.000. There is no separated 
region over the pressure surface for Re = 25,000. The location of 
separation point on the suction surface matches very well with the 
experimental data [25] for Reynolds numbers of 50,000 or less. Bons 
et al. [1] give no indication whether they observed any separated 
region over the pressure surface. 

We may note that the boundary layer over the entire blade is 
laminar, as evidenced in Fig. 3. This figure shows the dimensionless 
turbulent viscosity, p„ within a narrow range from 0 to 1.0 in the flow 
field for Re = 50.000 and B = 0 and 1 . Noting that the dimensionless 
laminar viscosity, p p , is nearly unity (dimensionless local temperature 
is near unity), it is clear that the effective viscosity within the entire 
boundary layer over the blade surface is essentially laminar. We may 
note the presence of separation on the suction surface in Fig. 3(a), 
while with the VGJs turned on, the flow gets attached (cf. Fig. 3(b)). 
Similar results were obtained for Re of 25,000 and 100.000. 

Figure 4 shows the velocity vectors, colored by static pressure, in 
the flow field around the blade for inlet Reynolds number of 50.000, 
and for B = 0 and 1. The thick separated region, with reverse flow, is 
clearly visible in Fig. 4(a) for the uncontrolled blade, while there is 
almost no separation on the suction surface for the controlled blade 
(Fig. 4(b)) with B = 1. On the pressure surface, the velocity profiles 
appear to be close to separation over a very small region, as discussed 
above regarding the results in Fig. 2. Similar results were obtained for 
Re of 25,000 and 100,000. 

Figure 5 compares the present computations with the experimental 
data [1] for the pressure coefficient on the baseline blade at two inlet 


Reynolds numbers and Tu = 1%. The comparison is very good on the 
pressure surface but not so good on the suction surface due to the 
presence of separation. While the experimental data exhibit the classic 
plateau indicating non-reattaching laminar separation, the analysis 
predicts significant diffusion through the separation region. Similar 
discrepancy was observed by Bons et al. [1] between their 
experimental data and the 2-D Navier-Stokes Vane Blade Interaction 
(VBI) code. 

Figure 6 compares the computed boundary layer profiles over the 
baseline (no VGJ) blade suction surface at three chordwise locations 
with the experimental data from Bons et al. [1] for Re = 50,000 and 
Tu = 1%. At 67% axial chord, the experimental profile was noted by 
Bons et al. [1] to be unsteady in the region y < 3mm. That is why the 
experimental data fall into a wide band in the near-wall region. At 
73% and 79% axial chord, the flow is separated, with the computations 
clearly showing reverse flow. For the experimental data, reverse flow 
direction could not be resolved [1] since a single hot-film probe was 
used. However, separate smoke visualization did confirm that flow 
reversal was present inside the separation bubble at the 73% axial 
chord measurement station. The comparison between the computed 
and experimental data is not very good, specially at 79% axial chord 
location, where the experimental data show much thicker boundary 
layer compared to the computations. This implies that the separated 
region is not well resolved by the computations. 

A similar comparison is shown in Fig. 7 for inlet Reynolds 
number = 100,000 and B = 0. Clearly at this high Reynolds number, 
the boundary layer is thinner compared to that in Fig. 6. However, the 
experimental data show a thinner boundary layer compared to that 
from the computations at all three chordwise locations. While the 
computations show separated (reverse) flow at both 73% and 79% 
axial chord locations, Bons et al. [1] report that the flow was close to 
separation at the 79% axial chord location but attached at the 73% 
axial chord location. Figure 8 compares the computed velocity profiles 
in the wake at 0.62 axial chord downstream from the blade trailing 
edge for two inlet Reynolds numbers against the experimental data of 
Bons et al. [1] with no VGJ. While the computations show little effect 
of Re on the wake width and relatively small effect on the wake depth, 
the experimental data show a larger effect of Re on both the wake 
width and depth. At Re = 50.000, the computed and experimental 
wake widths match but the computed wake depth is only about half of 
the experimental value. At Re = 100.000, the computed wake depth 
compares well with the experimental value but the computed wake 
width is about twice that of the experimental value. Such problems 
with computed wakes are well known, and point to the need for more 
work in this area. Moreover, in the present case, these problems are 
augmented by the large separation region over the suction surface of 
the blade. Any unsteady effects owing to the flow separation and 
wake shedding cannot be resolved by the steady code that was used for 
the present computations. 

Figure 9 shows a comparison between the presently computed loss 
coefficient, y, and the experimental values provided by Bons et al. 
[10]. While the experimental data show a sudden drop in loss 
coefficient around Re = 50.000, the computations show only a gradual 
decrease in the loss coefficient from Re = 25,000 to 100,000. This is 
again a reflection of the poor wake prediction noted above. Figure 10 
compares the computed boundary layer profiles over the blade suction 
surface at three chordwise locations with the experimental data from 
Bons et al. [1] for Re = 50,000 and B = 2. Both the computations and 
experiment show no separation (reverse flow). However, the 
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experimental data show steeper profiles than the computed results at 
all the stations, similar to the results in Fig. 7 for Re = 100,000 and no 
VGJ. Bons et al. [1] attribute the doubly-inflected (but attached) 
velocity profile at 79% axial chord station to the ultimate penetration 
depth of the vortex jets into the freestream. The computed results do 
not show any such anomaly. 

Figure 1 1 compares the computed velocity profiles in the wake at 

0. 62. axial chord downstream from the blade trailing edge for an inlet 
Reynolds number of 50,000 against the experimental data of Bons et 
al. [1] with and without VGJs. We may note that the data 
corresponding to B = 0 are also available in Fig. 8; they are included 
here to compare with the data for B 0. The computed wake with 
VGJs on (B ^ 0) is shallower and somewhat narrower than that without 
any VGJ (B = 0), as one would expect since the separation region 
vanishes when VGJs are on. However, the effect is not so dramatic, 
specially for the wake width, as for the experimental data. 
Experimentally the wake is very narrow when the VGJs are on 
compared to the case when they are not. The wake depth recovery, 
ratio of the wake depth for B ^ 0 to that for B = 0, is about the same 
for both the computational and experimental data. Also, there is little 
difference between the profiles corresponding to B = 1 and B = 2 for 
both the computational and experimental data. 

CONCLUSIONS 

The Glenn-HT code was used to compute the flow over a low- 
pressure turbine blade at very low Reynolds numbers with and without 
the use of vortex generator jets. The results are compared with the 
experimental data of Bons et al. [1] for steady VGJs. It is found that 
the code determines the proper location of the separation point on the 
suction surface of the baseline blade (without any VGJ) for Reynolds 
numbers of 50,000 or less. Also, the code finds that the separated 
region on the suction surface of the blade vanishes with the use of 
VGJs. However, the separated region and the wake characteristics are 
not well predicted. The wake width is generally over-predicted while 
the wake depth is under-predicted. Thus, there is a need to improve 
the prediction of wake characteristics, especially in situations where 
large separated regions may exist upstream of the wake, as in the 
present case. For resolution of the unsteady effects due to flow 
separation and wake shedding, use of an unsteady code is 
recommended. 
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Spanwise section of the multi-block 
viscous grid around the blade, with 
detail of grid over the blade span near 
and within the hole exits 
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Fig. 2 Computed shear stress on the blade surface at various inlet 
Reynolds numbers (Tu = 1%; B = 0) 




Fig. 3 Dimensionless turbulent viscosity within 0 and 1 in the flow field for 
Re = 50,000; Tu = 1%. (a) B = 0; (b) B = 1. 



Fig. 4 Computed velocity vectors, colored by static pressure, around the blade for Re = 50,000; Tu = 1%. 


(a) B = 0; (b) B = 1. 



Fig. 5 Comparison of pressure coefficient for 
baseline blade at two inlet Reynolds 
numbers and Tu = 1% 



u/mid-channel velocity 

Fig. 6 Comparison of boundary layer profiles of streamline velocity 

normalized by mid-channel velocity at three chordwise stations: 
0.67, 0.73 and 0.79 axial chord on the suction side. 

Re = 50,000; Tu = 1%;B = 0 



u/mid-channel velocity 

Fig. 7 Comparison of boundary layer profiles of streamline velocity 

normalized by mid-channel velocity at three chordwise stations: 
0.67, 0.73 and 0.79 axial chord on the suction side. 

Re = 100,000; Tu = 1%; B = 0 



Fig. 9 Comparison of loss coefficient vs. inlet Reynolds number for Tu = 1% 



u/mid-channel velocity 

Fig. 10 Comparison of boundary layer profiles of streamline velocity 

normalized by mid-channel velocity at three chordwise stations: 
0.67, 0.73 and 0.79 axial chord on the suction side. 

Re = 50,000; Tu = 1%; B = 2 
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y/blade spacing 

Fig. 8 Comparison of wake velocity profile at 0.62 axial chord downstream 
from the blade trailing edge for two inlet Reynolds numbers. Tu = 1% 
and B = 0 


y/blade spacing 

Fig. 11 Comparison of wake velocity profile at 0.62 axial chord downstream 
from the blade trailing edge for Re = 50,000; Tu = 1%; B = 0, 1, 2 
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